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Summary. — This paper offers a conceptually straightforward method for the cal- 
culation of stresses in polarisable media based on the notion of a drive form and 
its property of being closed in spacetimes with symmetry. After an outline of the 
notation required to exploit the powerful exterior calculus of differential forms, a dis- 
cussion of the relation between Killing isometries and conservation laws for smooth 
and distributional drive forms is given. Instantaneous forces on isolated spacetime 
domains and regions with interfaces are defined, based on manifestly covariant equa- 
tions of motion. The remaining sections apply these notions to media that sustain 
electromagnetic stresses, with emphasis on homogeneous magnetoelectric material. 
An explicit calculation of the average pressure exerted by a monochromatic wave 
normally incident on a homogeneous, magnetoelectric slab in vacuo is presented and 
the concluding section summarizes how this pressure depends on the parameters in 
the magnetoelectric tensors for the medium. 

PACS 02.40.Hw - Classical differential geometry. 

PACS 03.50.De - Classical electromagnetism. 

PACS 41.20.-q - Applied classical electromagnetism. 

PACS 41.20.Jb - Electromagnetic wave propagation. 

PACS 46.05.-f b - General theory of continuum mechanics of solids. 



1. Introduction 

The calculation of stresses in material media has extensive application in modern 
science. The balance laws of continuum mechanics offer an established framework for 
such calculations for matter subject to a wide class of constitutive properties that at- 
tempt to accommodate interaction with the environment in terms of phenomenological 
relations [1], [2]. Such relations are not always easily accessible via experiment, since the 
response of matter to internal and external interactions can be very complex. If one for- 
mulates these interactions in the language of forces derived from stress-energy-momentum 
tensors, then it is sometimes non-trivial to determine experimentally an appropriate 
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tensor that can be associated with a particular class of interactions on a macroscopic 
scale [3]— [8]. This problem has led to numerous debates over the last century about 
how best to formulate the transmission of electromagnetic forces in polarisablc media. 
Since the electromagnetic interaction is fundamentally relativistic in nature, the problem 
is compounded if one insists on a relativistically (covariant) theoretical formulation to 
compare with experiment in the laboratory. Judged by the large literature on this sub- 
ject, there is no universal consensus on how best to calculate forces in polarisable media 
and hence the needed experimental input into the subject has been of uncertain value in 
the past. However, modern technology — with the refined experimental procedures now 
available — offers the possibility that the appropriate constitutive relations for certain 
classes of polarisable matter can be determined experimentally [9] over a broad range 
of field intensities, frequencies and geometric configurations. Furthermore, new mate- 
rials with novel constitutive properties are being fabricated [10] and their response to 
time-varying electromagnetic fields also offers new potential for technological advances. 
With these points in mind, this paper offers a conceptually straightforward method for 
the calculation of stresses in polarisable media, based on the notion of a drive form and 
its property of being closed in spacetimes with symmetry. Section 2 outlines the nota- 
tion required to exploit the powerful exterior calculus of differential forms that is used 
throughout the article. Sections 3 and 4 relate the isometrics to conservation laws for 
smooth and distributional drive forms. Sections 5 and 6 discuss equations of motion in 
spacetime and how they may be used to define instantaneous forces on isolated domains, 
while section 7 deals with forces on domains with interfaces. The remaining sections 
apply these notions to media that sustain electromagnetic stresses, with emphasis in 
section 11 on homogeneous, magnetoelectric material. In sections 12—14, an explicit cal- 
culation of the average pressure exerted by a monochromatic, electromagnetic wave on a 
homogeneous, magnetoelectric slab in vacuo is presented and the discussion in section 15 
summarises how this pressure depends on the parameters of the magnetoelectric tensors 
for the medium. 



2. - Notation 

The formulation below exploits the geometric; language of exterior differential forms 
and vector fields on a manifold M [11]. Such a language is ideally suited to accommodate 
local changes of coordinates that can be used to simplify the description of boundary 
value problems and naturally encapsulates intrinsic global properties of domains with 
different physical properties. It is also makes available the powerful exterior calculus 
that facilitates the integration of forms over domains described as the images of chain 
maps and permits a clear formulation of notions such as energy, momentum, angular 
momentum, force and torque, by fundamentally relating them to isometrics of spacetime. 
In this framework, a p-form a belongs to SA^M, the space of sections of the bundle of 
exterior p- forms over M, while vector fields X belong to STM, the space of sections of 
the tangent bmidle over M. On a manifold with metric tensor g, we denote g{X, — ) by 

X G SA^M and conversely set X = X ior all X. In the following, a notational distinction 
between smooth (C°°) forms on some regular domain and those with possible singularities 
or discontinuities is useful. Smooth forms with compact support on spacetime will be 
referred to as test forms [12] and distinguished below by a superposed hat. Manifolds with 
dimension n will be assumed orientable and endowed with a preferred n-form induced 
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from the metric tensor field g. One then has [f 1] the linear Hodge operator * that maps 
p-forms to (n — p)-forms on M . If g has signature tg, one may write 

n 

(2.1) 3 = ^e'®e^r?y, 

where riij = diag{±\, ±1, . . . ± 1) and 

(2.2) ★l = eiAe2A...Ae", 

with tg = det{rjij) and {e*} a set of basis 1-forms in SK^M. The natural dual basis 
{Xi] is defined so that e^{Xj) = (5* and the contraction operator with respect to X is 
denoted ix ■ Covariant differentiation is performed with respect to the metric compatible 
Levi-Civita connection V, whilst Lie differentiation is denoted C 



3. Isometries and Drive Forms 

The notion of a drive form arises from the theory of gravitation in spacetimes M with 
isometries. In Einstein's theory of gravitation, the metric g of spacetime is determined 
by the tensor field equation 

Ein — T, 

where Ein e ST^M denotes the degree 2 symmetric divergence-free Einstein tensor field. 
Hence T must be a symmetric divergence-free degree 2 tensor field: 

V - r ^ 0. 

The tensor T is regarded as a source of gravitational curvature ! f^)! . If X is a Killing 
vector field generating a spacetime symmetry and * is the Lorentzian Hodge operator 
associated with g, then by definition 

Cxg^O 

and it follows that the drive 3-form 

TK = MT{K, -)) 

is closed on some domain Ij of M: 

dTK = 0. 



(^) The tensor T has dimensions of [MLT ^] (force) constructed from the SI dimensions 
[M], [L], [T], [Q] where [Q] has the unit of the Coulomb in the MKS system. 
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If the spacetime admits a set of Killing vector fields {Ki g ST Ij}, one has a conser- 
vation law for each Ki in every regular spacetime domain Ij [11], [13]. These may be 
supplemented with (tensor or spinor) field equations 

for all piecewise smooth (tensor or spinor) fields that interact with each other and 
gravity. These field equations may induce compatibility conditions and further (non- 
Killing) conservation laws 

(e.g. electric charge-current conservation). In phenomenological models, some of the 
field equations may be replaced by fixed background fields and source currents, together 
with consistent constitutive relations between these fields and currents. 

An observer field is associated with an arbitrary unit future-pointing timelike 4- 
velocity vector field U G STM. The field U may be used to describe an observer frame 
on spacetime and its integral curves model idealized observers. The drive form tk associ- 
ated with any K admits a unique orthogonal decomposition with respect to any observer 
frame U : 

where the spatial forms S SA^AI and € SA^M satisfy ijjp^ = ^uJk — ^- ^^3. 
local region, the conservation law dr^ = implies, in terms of the if-current j]^, the 
continuity relation in the frame U: 

If if is a spacelike tran slati onal Killing vector field and U a unit time- like (future-pointing) 
4- vector observer field 1(^1 then 

jU _ 

•JR = ~ t-utk 

is the linear momentum current (stress) 2-form in the frame U and 

is the associated linear momentum density 3- form in the frame U . If if is a spacelike 
rotational Killing vector field generating S'0(3) group isometrics, then j]^ is an angular- 
momentum current (torque stress) 2-form and is the associated angular-momentum 
density 3-form in the frame U . If if is a timelike translational Killing vector field, then 
is an energy current (power) 2-form and is the associated energy density 3-form 
in the frame U . In the following, attention is restricted to translational spacelike Killing 
vectors of flat spacetime and the computation of integrals of for a particular con- 
tribution to Tk associated with electromagnetic fields in homogeneous but anisotropic 



(2) The frame is inertia! if Vf/ = 0. 
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media of a particular kind. It will be argued that this formulation leads to a natural 
definition of integrated static forces in media with discontinuous material behavior and 
highlights the need for care in giving a practical definition of integrated force in media 
in the presence of time varying fields. 



4. Distributional Drive Forms 

To accommodate media with singular time-dependent sources of stress (e.g at surface 
interfaces or lines in space), introduce the distributional Killing 3-form tk^ on spacetime 
and its distributional source K,k^ satisfying the distributional equation 

(4.1) dTK^[/3] =/Ck''[/3], 

for all test 4-forms (3 [12] on spacetime. Consider a compact medium at time i, with 
spatial volume determined by the image of the spacelike i-parameterised immersion S'^t : 
Wz C ^ -M, evolving for a finite interval of time. Denote its immersed history in 
spacetime by the region Ii . Let I2 be a compact region of spacetime outside this medium 
history. It follows from (|4.ip that if tk^^ is the regular drive form in region Ji and tk^'^ 
is the regular drive form in region /2, then 

(4.2) dTK^^ =0 in /i, 

(4.3) dTK^^ =0 in /2 

(4.4) and ^^/(tk^^ - + JCk) = 0, 

at an evolving interface defined by the timelike, t-parameterised immersion T/^ s ■ S2 C 
M.^ —> M between Ii and I2 with a smooth interface drive form K-k on its image. The 
history of these images in spacetime is indicated schematically in figure (1). 

5. Equation of Motion for a Smooth Domain 

The notion of force (and torque) is implicit in the balance laws of classical Newtonian 
continuum mechanics. In the presence of time-varying fields, it is natural to associate 
energy, momentum and angular momentum with such fields in order to maintain the con- 
servation of these quantities for closed systems. The only sensible approach to defining 
force (and torque) density in such circumstances, where the balance law arises from the 
divergence of a total drive form for the system, is with respect to a particular splitting 
of this divergence. For systems without mechanical constraint, one assigns a smooth 
4-velocity V (and angular velocity) field to each smooth domain to describe the mo- 
tion. The jumps in these fields at interfaces between domains must be computed from 
(|4. 21 - 14. 4|l above. The 4-acceleration field A of each domain (and possibly its rate of 
change) will appear in one or more components of the split and the remaining terms in 
the divergence are often identified with total force (or torque) densities for the domain. 
However, unless one prescribes how to practically identify component contributions to 
the total force (for example by cancelling some of them by externally applied mechanical 
constraints) , there is no natural way to identify a canonical split of the divergence of the 
total drive form. In those situations where the interaction of matter and fields is station- 
ary or static, one can appeal to static experiments with non-moving media to try and 
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Fig. 1. - The partition of spacetime M by the history of a compact medium (with boundary 
E^s U E^jjj U E^t), evolving with 4-velocity V . The timelike vector field U defines a frame, N is 
a unit, spacelike vector field and K \s a, Killing vector field. 



give an unambiguous definition to material body forces. For electromagnetic interactions 
with polarisable media, comparison with experiment is difficult, since the choice of drive 
form is very model dependent for many materials. However modern technology — with 
the refined experimental procedures now available — offers the possibility that the appro- 
priate constitutive relation for certain types of matter can be determined experimentally 
over an extended parameter range. [10] . 

To illustrate these general remarks, consider an uncharged (unbounded) medium con- 
taining a fixed number of constituents, with number density A'' G SK^M and^mass density 
p = moAf, mo > 0, in Minkowski spacetime with mass conservation d{p-kV) = 0. Write 

^ V . ^ field 

where V is the (future-pointing) unit, time- like 4-velocity field of the medium and K a 
Killing vector field. For a simple medium with a smooth mass density, suppose 

Tx^= co''pg{V,K)^V, 

with Co the speed of light in vacuo, then 



dTK = 
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yields the local equation of motion [11] for the field V: 

where fx = -kdrR^''''-'^ and the 4-acceleration 1-form A = VyF. If VF = 0, then A = 
and the motion of the medium is geodesic. The medium is then static in the frame where 
U = V. Contracting the local equation of motion 

Co V^W^l + t^Tif-^*^" = 

with the observer field U and integrating over the volume S^t yields 

K t\ — JkI^ t\, 

where 

J 

Here the mass 3-form 

;U[/ = - Co ★ [/ 

and the total instantaneous integrated K-dnve component on E^t at time t in the U frame 
is 

/^[S^]^ / iudTK^'^^'^. 
J 



6. — The General Integrated Force Form on a Regulcir Domain Ij c M 

In Minkowski spacetime, one has a global basis of parallel unit space-like transla- 
tional Killing vector fields {Ki, K2, K^). In local Cartesian coordinates {t,x,y,z}, with 
g = — Co ^dt ® dt + dx ® dx + dy ® dy + dz ® dz: 

i^. = ^,i^. = ^,i^3 = f 

ax ay oz 

One can then define the instantaneous integrated force 1-form on S^t at time t in the U 
frame to be 

(6.1) f[i:\]^j2fK,[T.\]%. 
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Then, if N is any unit space- like vector field on E'^t, the instantaneous integrated force 
component in the direction TV acting on S'^t is 

(6.2) r[I]^](iV)^^/^^.[I]^]i?;(7V). 

In an arbitrary (possibly non-inertial) frame U and domain Ij C M 

where 

(6.3) f'"^'^ = -iuTK^'"^"^ 

is the total Cauchy stress 2-forn Jf^ . If one identifies an electromagnetic if-drive tk^^^^^ 
in Tx-''*''" , such that 

and 

„ U field /,■ „ U EMI, , „ U rem I, 

one then has 

prnech UI, p3^] ^ pEM UI, p3^] + premUI, p3^] ^ fEM UI, p3^] ^ jrem UI, p3^] ^ 

where 

Tjrem U In rv3 1 — / r ^ rem J. 



pEMUI,^^3^^^_ f jr^^^ 



EMI, 



denote integrated rates of change associated with field momenta in tk^'^"^ and tk^^^^^ 
respectively, and 



U rem Ij 



denote integrated forces associated with stresses in tk^"^"^ ^' and tk^^^^ respectively. 



('') This follows from the Cartan identity: Cx — ixd + dix for any X G STM. 
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7. — The General Integrated Force Form in an Irregular Static Domain Com- 
posed of Different Media 

Suppose E^t = Ij with the Cauchy stress 2-form for domain Ij in 

a Minkowski spacetime with frame U = In general, tk^^^'''^ must contain (time 

dependent) constraining forces to maintain the overall equilibrium condition 

7 field r\ 

I 8 dri. =0 

from stresses in each sub-domain Ij of E^j. In the (possibly constrained) static case, 
^utk = and each 2-form 

contributes an integrated reaction force on E^t from domain Ij . 

In general, each 4-vclocity V^^ G STIj must be determined from the jump conditions 
for tk^^^^'^ ■ In the static case, one has all = with V^^ = U, and one may define 

the net integrated K-foice for E^t in the frame U: 

[E^] = y [ i^dril''" =y I da'i ^'''"'^ =y [ a'i ^'''"'^ . 

There may be additional sources of stress with support on submanifolds of M. Singular 
sources of stress in the electromagnetic field include charges, currents and their multi- 
poles, with support on points, lines or surfaces in space [7]. If the integrals on 
the right below are finite, the most general integrated force can then be written so as to 
include such distributional sources: 

/^PM ^ E / 4'""'^ +E / 4'"'''^ +E / 7k 

in terms of line stress 1-forms jk and surface stress 2-forms kk- 

A number of sources of interfacial stress depend on the local mean curvature normal 
of the interface. For example, if the history of the interface d Ij is the spacetime hyper- 
surface / = with unit spacelike normal A'' = then the scalar (Tr H) is defined 
by 

diNiu * 1 = {Tr H) iu*l 

and 77 = (Tr H) N is the mean curvature normal. Surface tension at an arbitrary in- 
terface depends on r} and the local surface tension scalar field 7, yielding the particular 
interface forces: 

f U fieldii f . . . T 

/ 7k = / '^iNtUlK 
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8. Electromagnetic Fields in Spacetime 

Maxwell's equations for an electromagnetic field in an arbitrary medium can be writ- 
ten 

(8.1) dF = and d * G = j, 

where F e SA'^M is the Maxwell 2-form, G £ S A^AI is the excitation 2-form and 
j G SA^M is the 3- form electric current sourc jf^ . To close this system, "electromag- 
netic constitutive relations" relating G and j to F are necessary. The functional tensor 
relations 

G = Z[F] 

and 

are typical for idealized material without electrostriction losses. 

The electric 4-current j describes both (mobile) electric charge and effective (Ohmic) 
currents in a conducting medium. The electric field e S SA^AI and magnetic induction 
field b G SA^M associated with F are defined with respect to an observer field U by 

(8.2) e — iuF and Cfjb = iu -k F. 
Thus, if/e — ijjb — and with g(U, U) = —1, 

(8.3) F ^ eAU -i^icabAU). 

Likewise the displacement field d £ SA^M and the magnetic field h £ SA^AI associated 
with G are defined with respect to U by 

(8.4) d^iuG and —^iu*G. 

Co 

Thus, 

(8.5) G = dAC7-* AC/^ , 



(*) All electromagnetic tensors in this article have dimensions constructed from the SI dimen- 
sions [M], [L], [T], [Q] where [Q] has the unit of the Coulomb in the MKS system. We adopt 
[g] = [i^], [G] = [j] — [Q], [F] — [Q]/[eo] where the permittivity of free space eo has the di- 
mensions IQ^ T^M~^ L~^\ and co = ^^^^^ denotes the speed of light in vacuo. Note that, with 

[g\ — [F^], for p-forms a in n dimensions one has [-ka] — [a][L"^^'']. 
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with iijd = ijjh — 0. The spatial 1-forms e, 6, d, h are fields on a general spacetime 
defined with respect to the frame U, which may be non-inertial if dU ^ 0. 



9. Time-Dependent Msixwell Systems in Space 

In the following, attention is restricted to fields on Minkowski spacetime. This can 
be globally foliated by 3— dimensional spacelike hyperplanes. The Minkowski metric 
on spacetime induces a metric with Euclidean signature on each spacetime hyperplane. 
Furthermore, each hyperplane contains events that are deemed simultaneous with respect 
to a clock attached to any integral curve of a future-pointing, unit, time- like vector field 
U = defining an inertial observer on Minkowski spacetime and the spacetime Hodge 
map ★ induces a Euclidean Hodge map # on each hyperplane by the relation 

*1 = Co A #1 = #1 A U. 

The spacetime Maxwell system can now be reduced to a family of parameterised exterior 
systems on R'^. Each member is an exterior system involving forms on depending 
parametrically on some time coordinate t associated with U. Let the (3 + 1) split of the 
4-current 3-form with respect to the foliation be 

j = -JAdt + p#l, 

with i a J = 0. Then, from ((8?T|) . 

(9.1) dj = 
yields 

(9.2) dJ + /3#l = 0. 

Here, and below, an over-dot denotes (Lie) differentiation with respect to the parameter 

t (d = £ a a for all a) and d denotes exterior differentiation on such that 

at 

d = d + dt A Ca_. 

at 

It is convenient to introduce on each spacetime hyperplane the (Euclidean Hodge) dual 
forms: 



E = #e, D = #d 

B = #b, H= #h, J = # J, 
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SO that the (3 + 1) spUt of the spacetime covariant Maxwell equations (|8.ip with respect 
to U = —Codt becomes 

(9.3) de = -B, 

(9.4) dB = 0, 

(9.5) dh = J + i?, 

(9.6) dD = p#l. 

All p-forms {p > 0) in these equations are independent of dt, but have components that 
may depend parametrically on t. 



10. Electromagnetic Constitutive Tensors for Linear Media 

Attention will now be turned to integrated electromagnetic forces on a class of polar- 
isable media. This requires a discussion of a class of electromagnetic constitutive tensors 
for linear media. In general, the excitation tensor G is a functional of the Maxwell field 
tensor F and properties of the medium 

G^Z[F,...]. 

Such a functional induces, in general, non-linear and non-local relations between d, h 
and e, b. Electrostriction and magnetostriction arise from the dependence of Z on the 
elastic deformation tensor of the medium. For general linear continua, one may define a 
collection of constitutive tensor fields Z ^''^ on spacetime by the relation 

G = Sf^oZW[VF,...], 

in terms of the spacetime connection (covariant derivative) V. 

In idealized (non-dispersive) simple media, one adopts the simplified local relation 

G = Z{F), 

for some degree 4 constitutive tensor field Z and in the vacuum G = cqF. Regular 
linear isotropic media are described by a bulk 4-velocity field V, a relative permittivity 
scalar field e,. and a non- vanishing relative permeability scalar field /i,- . In this case, the 
structure of Z follows from 

G ~ 

— =er ivF AV - fi:;^i^(iv^F AV) = ivF AV + —F. 

In a comoving frame with U ~ V , this becomes 

d = eoer e and h = {^ol^rj^^b. 

To discuss linear (non-dispersive, lossless), inhomogeneous, anisotropic media, it is con- 
venient to describe Z in a particular basis associated with the medium. Since Z is a 
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tensor that maps 2-forins to 2-forins, in any spacetime local frame {e^ ,e^}, one 
may write 

iGabe" A = \Z''\^F,ae'' A e^ 

where 

^ ab — ~^ ba — ~^ ab — ^ ba- 

Thus, Z can be described in terms of spatial rank 3 tensors on spacetime, relating 
observed electric and magnetic fields in some frame U, with: 

h = C^<^(e) + C^b(6). 

In such a frame, the medium is said to exhibit magneto- electric properties in general. 
If and are non-zero in the co-moving frame of the medium, it is called magne- 
toelectric. If and (^^^ are zero in the co-moving frame of the medium, it is called 
non-magnetoelectric. The spatial tensors and C,^'' may be non-zero in a non-comoving 
frame for a non-magnctoolcctric incidium. Duo to the behaviour of electric and magnetic 
fields under Lorentz transformations, all materials exhibit magnetoelectric properties 
in some frame. Thermodynamic and time symmetry conditions impose the relation 
Z = Zt [14] or 

^det ^ ^do ^ ^hbt ^ ^hb ^dbt ^ _^he 

in all spacetime frames, where the adjoint T'^ of a tensor T which maps j3-forms to p-forms 
is defined by: 

a A *T(/3) = /? A ★T^ (a) for all a, /3 G SA^M. 



11. — Homogeneous Dispersive Magnetoelectric Media 

In dispersive media, constitutive relations between the spatial fields e, b, d, h are non- 
local in spacetime. If the medium is spatially homogenous, so that it has no preferred 
spatial origin, then it is possible to Fourier transform the fields with respect to space and 
time, and work with transformed local constitutive relations. 

For any real valued p-form a, define its complex valued Fourier transform dk,ui by 

/OO POO 
du) dk cek,oj expi{k ■ i — ut), 
-OO J — OO 

where fe e R^. Then the source free Maxwell system reduces to 

(11-2) KAe^^^=u;B^^^ 
(11-3) KAhk,^ = -u:Dk^^, 
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where the real propagation wave l-form K — k dr E SA^M. The remaining transformed 
Maxwell equations KhB^.^^ = and KhD^^^ = follow trivially from (fTT^ and pT3)) 
when w 7^ 0. It also follows trivially that uj^^k = ( i-^- ^k w perpendicular to ^ 
Similarly, B^^ A K = and D^^ AK = 0. 

We assume that the magnetoelectric constitutive relations take the form 

(11-4) dk,u; = CtA^k^uj) + 'ctA^k,J 

(11-5) = Cfc^c^(efe,c^) + C^A^k.J- 

These will (by convolution) give rise to non-local spacetime constitutive relations. We also 

, , . , , , rde -rdh rhe 

mamtam the above symmetry properties on the magnetoelectric tensors Qk uj^^k un>.k un>.k 
Substituting \11A\ and (|11.5|) in (|11.2[) and (jll.Sp yields a degenerate l-form linear eigen- 
equation for e^, ^: 

k.uj (^fc,cj ) j 

(11-6) +#(KACfc!'.(#(fsrAefc,J)) =0. 

The field bf,^ then follows from (|11.2p . (up to a scahng) and dj^^, hf,^ from (|11.4p . (|11.5p 
respectively. Equation (I11.6P may be written 

(11-7) -DkAekJ^Q, 

defining the 1 — 1 tensor T>k.Lu- For non-trivial solutions e,^^, the determinant of the 
matrix Vk.u representing 'Dk^cj must vanish: 

(11.8) detiVk,^) = 0. 

Note that, in general, the roots of this dispersion relation are not invariant under the 
transformation K ~K. If one writes fc — k.\k\ in terms of the Euclidean norm |fc|, 
and introduces the refractive index M = 1^1"^ > ^ind k in place of fc, then solutions 
propagating in the direction described by fc with angular frequency > correspond 
to roots of pT8)) (labelled r) that may be expressed in the form Mr ~ J-'r{k,Lu). Thus, 
there can be a set of distinct characteristic waves each with its unique refractive index 
that depends on the propagation direction fc and frequency lu. When the characteristic 
equation (Ill.Sp is a quadratic polynomial in A/"^ and has two distinct roots that describe 
two distinct propagating modes for a given w, the medium is termed birefringent. Roots 
Af^ such that Afr{k,uj) ^ A/'r(— fc,w) imply that harmonic plane waves propagating in 
the opposite directions ±fc have different wave speeds. 

Each eigen-wave will have a uniquely defined polarisation obtained by solving the 
independent equations in ()11.7p for ^, up to normalisation. Since ^ is complex, it 
is convenient to introduce the eigen-wave normalisation by writing 
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in terms of the complex 0-form ^ and complex polarisation 1-form fi^ ^ , normalised 
to satisfy 

(11-9) ^A#n^,^ = #l 

for each r. If one applies # e£ ^ A # to (|11.6p . making use of the symmetries between the 

rde ^ dh -he r^hb . . , , . 

real magnetoelectric tensors Ck.ui Ck.ujj Ck.uji Cfc,i^i and evaluates it with the eigen-wave 
(J I one obtains the real 0-form dispersion relation for the characteristic mode r: 

+^ # A K A Cfc^JnL, j) + # A K A Cfc'. A nlj)) = 0, 

where K — —Afk ■ dr in terms of M and k. 



12. Electromagnetic Stress-Energy-Momentum Tensors 

There has been intense debate over many decades about the appropriate choice of 
electromagnetic stress-energy-momentum tensor that transmits forces in a (moving) po- 
larisable medium [15]. In 1909, Abraham introduced the symmetric electromagnetic 
stress-energy-momentum tensor T^^^ for a medium with 4-velocity V: 

2'j-EM ^ -i^F ®eG -iaG^i^-F -i.{F M.G)g + V ®s + s®V, 

where ia = «x„, = 9°^'^ ib in any vector basis {Xa} and 

1 



s = ★ ( -^e^ A /i^ A y - Co A A F ) , 



where ^ iyF, cq 6^ — iv * F, dX — lyG and = iv*G, 



Co 



are fields defined relative to the motion of the medium, so that 



G^d"" AV-* (^^AV^ 

with 

G = Z{F). 

For any Killing field K the drive form associated with Abraham's electromagnetic stress- 
energy-momentum tensor is 

(12.1) t|*^ = i (^F AiK*G-iKGA-kF + s{K)*V + V{K)irs 
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It follows from MM. and KB that 



(12.2) = = i ( e W#d + d(i^)#e + h{K)#b + b{K)#h) 



-i# e A #d + 6 A ifhj #K + \u{K) ^-^e Ah+ cadhb 



-^iu (^K Aiv * sj - V{K)iu * s 



and 



(12.3) = --t^W (b A #/i + e A #d) + i f — e A + cq d A b ) A 

2 2 V Co y 

-]au [k hU Mv^sj +y(-ftr)i[/(*sA ?/), 

where iiT-L = K + U{K)U. 

By contrast, Minkowski (1908) introduced the non- symmetric electromagnetic stress- 
energy-momentum tensor T^^^ where 

(12.4) T^^ ^ -~iaF ® i'^G - ^ ★ (F A ★ G)(7, 

which exhibits no explicit dependence on the medium 4-velocity V. The corresponding 
drive form is 

and (|6.3[) . (|8.3p and (|8.5p yield in this case 

(12.5) EE cr|^ = h{K)H^b + e{K)ifd + —U{K)e A 

Co 

-i#(eA#d+bA#/i) 

and 

(12.6) = Co d A b A - ^U{K) {eh#d + bh#h). 

More recently other choices for an electromagnetic stress-energy-momentum tensor 
have been proposed which in themselves simply imply different constitutive relations [16] 
with respect to a particular total stress-energy-momentum tensor. In [17] , [18] , it has been 
argued that different choices of the electromagnetic stress-energy-momentum tensor for 
linear polarisable media are equivalent to different choices of Z and a different partition of 
the total stress-energy-momentum tensor for the computation of so called pondermotive 
forces that arise from the divergence of terms in its decomposition. Furthermore, it was 
shown how particular choices of the dependence of Z on the gravitational interaction led, 
via a covariant variational formulation, to either the Abraham tensor or a symmetrized 
version of that proposed by Minkowski. 
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In the following, we illustrate how the general theory of drive forms outlined above 
offers a natural tool to discuss the computation of particular electromagnetic forces for 
materials that exhibit magnetoelectric properties (at rest) in the laboratory, for a par- 
ticular choice of electromagnetic drive form. This is an essential step in any program 
that attempts to confront experimental measurements of such forces with theoretical 
prediction. 

To facilitate this calculation, an electromagnetic drive form associated with the tensor 
obtained by symmetrizing p2.4p will be chosen: 

(12.7) TK^'^'^ :^^{FAiK*G-iKGAirF). 

It follows from (16.31). (18.31) and (18.51) that with this drive- form: 



(12.8) Jj^ = = i ( e(i^)#d + d{K)#e + h{K)#b + b{K)#h 



( e A #d + 6 A + \u{K) (^-^e Ah + codAb 



and 



(12.9) Pk = 7; 1 — e Ah+codAb]AK^- lu(K) (e A #d + b A #/i) . 
2 V Co J 2 

For a medium at rest in the laboratory, U = V = -^dt- Furthermore, if U{K) = 0, 
the 2-forms (|12.2p and (|12.8p coincide, so the following analysis does not discriminate 
between the choice of tensors (|12.ip and (|12.7[) . However, in this case the instantaneous 
densities (|12.3I) and (|12.9p are different. But, for the polarised monochromatic plane 
waves discussed below, the time-averaged tensors based on (|12.6p and (|12.9p also coincide. 
If the fields are all differentiable in the medium described by p2.7p , one readily obtains 

drf*^ = ]^{iKdG A-*:F + iKG Ad-*:F - F AiKdi^G) 

= FA«Kd(|) - (f+ (^)) A ^KJ+GA^.d. 

where 



dF = 0, dirG = j, G^cqF + U, €„ d F = j - d i^U. 

Thus, non-zero bulk integrated static electromagnetic forces from such fields require 
dH 7^ 0, d n 7^ (magnetisation or electrical polarisation inhomogeneities) or j ^ 
(non-zero local source current or charge density). For a neutral homogeneous material 
therefore, we consider a medium whose electromagnetic properties change discontinuously 
at some interface. 
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13. The Magnetoelectric Slab 

In t erms of the rank 3 identity tensor Id in space, consider an infinitely extended 
slal jf^)! of magnetoelectric material with 

- de 

(13-1) Ck,iu = ^k.uj Id 

(13.2) C't = 

The slab has width L and parallel interfaces (with the vacuum) at a; = and x — L. It 
is oriented in the laboratory frame {^x^^y,^z\^ so that C'k^ takes the particular form 

(13.3) (fk,ui = Pi,k,uj dz®dy+ I32,k,uj dy®dz, 



Incident, E^j^ 




Reflected, E'^ 
< 



ly 



z 



z 
1 1 



Transmitted, E^J^ 



Reflected, E^^ '^ 
< 



Transmitted, E^^^ 



III 



x = {) 



X = L 



Fig. 2. - Geometry of the magnetoelectric slab and the electric field amplitudes in the three 
regions. 



In this frame, the modes associated with the branch ()13.6p of the dispersion relation 
below will be polarised in the direction and those associated with the branch (|13.7p 



db 



will be polarised in the direction dy . The matrix representing C,k u takes the form 



db 



(13.4) [Ck 




(^) Such a medium has been considered by Hehl and Obukhov in their classical analysis of the 
Feigel effect [19]. 
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It follows that 

- he 

(13.5) Ck.uj = -l32,k,u dz®dy- Pi,k,Lj dy (g) dz. 

With this choice of orientation of the slab, the spatial region < x < L will be denoted 
// and the region with x > L denoted ///. The electromagnetic fields induced in its 
interior by a plane monochromatic wave normally incident from the right (in region /, 
X < 0) propagating in the direction dx with polarisation in the direction dy can now be 
readily determined. 

From (jll.Sp . the dispersion relation associated with one polarised eigen-mode of e,^^^ 

is 

(13.6) efc,w/ifc,wW^ - fc^ - 2(3iM,uikuj = 0, 

while that associated with the other polarised eigen-mode of ^ is 

(13.7) efc,wA*fc,w^^ - fc^ + 2(32M,uikuj = 0. 

Each relation can describe propagating modes with angular frequency w > moving in 
a direction determined by sign{k) dx with phase speed l^^/fcl- Since this ratio depends 
on the values of Pi^k,uj or P2,k,ui, it may exceed the speed of light in vacuo. In principle, 
such modes can contribute to the synthesis of wave packets. However, in the following, 
we restrict to monochromatic incident waves and work with constitutive parameters that 
inhibit super-luminal waves, with real constants ek^u = e > 0, fik,u = /i > 0, /3i,fc,t^ = 
P2,k,uj = /32. For an incident wave with complex amplitude £, no loss of generality 
arises by taking lo > and writing the solution ^: 

(13.8) ^fc^ — ^ (^exp{ik^j^ X — icu t) dy + E^^^ exp(ik^j^x — iuj t) dy^ 

(13.9) ejj^ — £ (^E^J^ exp{ik^J^ X — iujt) dy + E\^^ eyij>{ik]^^ x — iuit) dy^ 

(13.10) e["j = EE^J^yeM^k'n'^x - i^t) dy 

where k^^^ denotes a real root of the dispersion relation (|13.7p associated with the po- 
larisation eigenvector dy with sign{k^Jy) > 1, describing a polarised right-moving wave 
in the slab (region II) . Similarly, fc^^^ denotes a real root of the dispersion relation as- 
sociated with the polarisation eigenvector dy with sign[k^j^y) < 1, describing a polarised 
left-moving wave in the slab (region //). In general, these wave numbers are different. 
In the vacuum regions, fcjj^ = —k^j^ = k^j^^^ = ui/ cq . 

If JIq (ilj^) denotes the pull-back of forms to the interface a: = (a: = L), the interface 
boundary conditions [20] are 

(13.11) {et ~ C') - (C' - C) = 

(13.12) (hii - hi'::) - ^1 {hk.: - hZ') = o 



(13.13) % (bH - Bi':) = ni (b'^: - Bi"j) = 

(13.14) n* (bil - Di'/) = ni \b'l: - bZ') - 0, 
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yielding the linear system for the dimensionless complex amplitudes E^j^^ , eU ^ , e\^ ^ , eU^ ^ : 



Jiy, 



i + E[y = E'^y + E['y 



Jiy 



(32 1 E'Jyexp{ik'j(yL) + 



Jiy 



jJLUJ 



P2]Ei'yeMikrL) = 



ki"yEU'y 



This system of equations has the solution: 



E{y = 



E 



II y 



(13.15) 



E 



II y _ 



E 



Illy _ 



n 


'■■;;■) ( 


II I'll 'h^' - 




(wpf//, 




M (^•;;' - 


'^■) ( 






exp(/7, 






- '1" 




k['A <.xp(/(/,-|/" 


+ A-f " 





^Jiiy 



where it is convenient to introduce 

= (exp(zfc^^^L) - exp(^fcf ^L) ) [^iV' (mo/32C. + ^^j"^) ± ^l^ll(i20J (/JaM^ - 4'' " ^"') 
+110 {j32UJii^k^R^y + ^o^fi ^^"^) ± ( fcfl ^ exp{ik[^'"L) - k"^ exp(iA;^^^L) j M/^o^j 
+ ^ k^j^y exp{ik^j^y L) — k^j^y exp{ik^^^y Lfj /i/io^j 

with 

k^+' = k^ii and 



4^ 



With the electric field amplitudes determined, the complete set of polarised fields 
{^ku}^^k,oj^^k,ui^^k,u>} e&ch region is determined. For completeness, these fields are 
given in the appendix. 



14. — Average Pressure on the Magnetoelectric Slab 

To calculate the average pressure on the sides of the magnetoelectric slab, one inte- 
grates the Maxwell- Cauchy stress 2-form over the 2-chain (surface) = 
indicated schematically in figure (3). The image of is the boundary of a box of height 
H, width W and length L, with faces JIq and fl^ in regions / and /// respectively, par- 
allel to the surfaces of the slab. Integrating over a box with faces wholly within II would 
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give zero total force, since region // is homogeneous. Since the fields are independent of 
z, contributions to the integral from the oriented chains fii and Q.2 cancel. 



I 



— 

f2i 



9.2 



// 



a; = 



X = L 



H 



I ± 



/// 



Fig. 3. - Geometry of the 2-chain Q, used to calculate the time-averaged integrated pressure on 
the magnetoelectric slab. 



The above fields yield a net pressure on II that fluctuates with time, with a non-zero 
average. If A{r, t) is a scalar field, its average over any time interval T is 

1 r 

<A>{r) = -j^ A(r,t)dt. 
Hence, \{ B{r, t) is another scalar field. 



<AB>{r) = l- I A{r,t)B{r,t)dt. 



Furthermore, if 

A = Re ( A{r) exp{-iwt)^ € SA^M 
B = Re ( B{r) exp(-ia;t)) e SA'^M, 

where A, B are complex, then 

AaB = ^Re(^AABexp{-2iujt)^ + ^Re(^AS), 

so 

(14.1) < AAB > (r) = ^Re(^AB), 
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if we take T = 2^. Thus, gives 
< 



Furthermore, with U = -^dt and K ~ d^, the a;-component of the time-averaged 
MaxweU-Cauchy stress 2-form is 



-i#Re (el^ A + bl^ A ^IQ 



k,L 



which reduces to 



(14.2) < ao^ > (r) = - J#Re [el^ A #dl^ + bl^ A #hlS) dyhdz. 

The time-averaged integrated force is given bj{{£}| 

Denote the time-averaged stress forms due to the fields in regions / and /// by 

<(t])^>^ a' dy A dz, a' = - J#Re (^e,^^ A #4;^^ + bil A ifhil^ 
< ag^ > = a"^ dy A dz, a'^' = -i#Re (e^^^^ A + ^ A 

Thus, the time-averaged net force on the magnetoelectric medium contained in the region 
bounded by is 

< fl^^M > = {n*a^ - nW") / / dydz = {%a^ ~ ^W") A, 

Jo Jo 

where A — WH and the time-average integrated pressure < [ri] > = — . 

(^) The minus sign occurs due to the opposite orientation of the opposite faces of fl. 
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Calculating the pull-backs of 



£0 1 + 2Re(i?[^exp(z[4^ ~ k[y]x)) + \E[y\' 



1 



Jy\2\ j^iy\2 



(a 



yields 



{i + \E[y\- 



since k 



ly 



■Jy 



iJiiy 



Since the time-averaged body force < Cup^ >= 
for given by (|12.3p . (|12.6p and (|12.9p . it follows that the average pressure on the 
magnetoelectric slab is given in terms of the solution (jlS.lSp by 



(14.3) 



< PxM > 



(I 



E 



iiiya 
R I 



15. Conclusion 

The magnitude and sign of < Pxi'^] > depends on e = erCo, P = PrPa, Pi and /32, where 
Co — , ^ As noted above, the wave numbers ki^^ , klJ^ that follow from the dispersion 
relation determine the nature of the propagating wave in region //. For the case under 
discussion here, where the parameters e^., Pr, f3i, P2 are constant, it is of interest to write 
the dispersion relations in terms of the dimensionless ratio of the wave speeds w = 
where v = j, vq = and the dimensionless parameters bi = —(ii/y'Tp, 62 = (32/y/ep'- 

+ 2hiw -1 = 
+ 262W -1 = 0. 

Then the sub- luminal condition |-^| < 1 implies \w\ < ^trPr- The relation between w 
and either h\ and 62 can then be seen from the relation of the two branches of the loci 
where the expression 'uP' -V 2bw — 1 vanishes in the w-b plane. For a; > 0, values of w 
in the upper (lower) half plane correspond to left (right) moving waves. Furthermore, 
propagating sub-luminal monochromatic waves will only occur in // for real b, yielding 
real values of w in the range —^/erPr < w < ^CrPr- It is clear from these considerations 
that the relative sign between /3i and (32 can have a significant effect on the behavior 
of the propagating modes in the region // and hence on the nature of the force on the 
magnetoelectric slab. 

The authors feel that the approach adopted in this paper for the calculation of static, 
time averaged and instantaneous forces, offers a conceptually unambiguous method of 
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considerable generality. Once one decides on the drive form appropriate for any subsys- 
tem in interaction with external fields, it has immediate application to moving media 
(in arbitrary relativistic or non-relativistic motion) and can be extended to matter with 
material losses. Work is in progress to extend the methodology to inhomogeneous media 
with more general constitutive properties and this will be reported elsewhere. 



* * * 
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Appendix 



Electromagnetic Fields in the Three Regions For a y-polarised harmonic electro- 
magnetic wave with angular frequency u; > 0, incident normally from the left on a 
fixed magnetoelectric slab, the electric field solutions in the three regions are given by 
(|13.8p — (|13.10p . For completeness, the remaining fields in these three regions are given 
here. The magnetic induction fields follow from (|I1.2p : 



-7 £ k^^ Ek^^E^^ 

^k'uj — ~ exp{ik^^ X — iuj t) dz H — — exp{ik^j^x — iuj t) dz 

^ UJ ' UJ 

-II Ek^^^E^^^ Sk^^^E^^^ 

bfc J — — exp{ik'[^^x — ioj t) dz H — exp{ik^]^^^x — ioj t) dz 

-III £ k^^^^E^^^^ 

^k,J^ = ~ ~ — exp{ikl^'^x — iujt) dz. 



The electric displacement 1-forms in regions / and /// are given by the vacuum consti- 
tutive relation dk uj — ^o^^^, whereas the electric displacement 1-form in region // is 

given by the constitutive relation (lll.4|) , with the spatial tensors C'kuj ^-^^ (f^uj given by 
(|13.ip and pS.Sp respectively: 



^feL — ^ exp{ik^j^x — iuit) dy + £ foE^j^^ exp{ik^j^x — iut t) dy 

^^k<!j = £ i £ + P'^^R I jT^i^ y exp(ifcjj^^a; — iujt)dy + £ ( e + ^"^^^ | e\^ ^ exp(ifc^^^x — iuj t) dy 
^k.uj^ — £ £oEj^'^ exp{ik^J^^x — iujt)dy. 



Similarly, the magnetic I-forms in the regions / and /// are given by the vacuum con- 
stitutive relation — f-^Q^bj^ ^, whereas in region //, the magnetoelectric constitutive 

relation pi .5p , with the spatial tensors c'k'^ and c'k^ given by (|I3.2p and (|13.5p respec- 
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tively yield 



fly 



hk.!!> - ^ [—- P2\E^jly expiik^J^x -iujt)dz + £ ( ^ - ) ^ exp(z/fcf - luo t) dz 



c Uy c ijy pi y 

exp[ikj^x — lU! t) dz H exp[ikjfx — ilo t) dz 



Mow 



Mow 



' l,iiy 

ft, n 



kV^ 



Illy j^IIIy 



h 



Illy _ £k^__rEj 

fJ-OLO 



{ik^J^^x — iu! t) dz. 
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